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ABSTRACT 

The P1 promoter of the galactose operon in 
Escherichia coli is one of the best studied examples 
of 'extended -10' promoters. Recognition of the P1 
promoter does not require specific contacts be- 
tween RNA polymerase and its poor -35 element. 
To investigate whether specific recognition of the 
-35 element would affect the regulation of P1 by 
GaIR, we mutagenized the -35 element of P1, 
isolated variants of the -35 element and studied 
the regulation of the mutant promoters by in vitro 
transcription assays and by mathematical modeling. 
The results show that the GalR-mediated DNA loop 
is less efficient in repressing P1 transcription when 
RNA polymerase binds to the -10 and -35 elements 
concomitantly. Our results suggest that promoters 
that lack specific -35 element recognition allow 
decoupling of local chromosome structure from 
transcription initiation. 



INTRODUCTION 

Most Escherichia coli promoters transcribed by the major 
RNA polymerase (RNAP) (a 2 (3P'a 70 ) contain two con- 
served sequence elements located at positions 10 and 
35 bp upstream of the transcription start point (tsp) 
(1,2). These promoter elements are recognized by the rj 0 
subunit of RNAP. The conserved regions 2.4 and 4.2 of 
a 70 recognize the —10 and —35 promoter elements, re- 
spectively (3,4). However, a minor class of promoters 
lacks a recognizable —35 element. These promoters are 
contacted with a different pattern, where region 2.4 
binds to the —10 hexamer and region 2.5 makes sequence 
specific contacts with a short sequence located 1 bp up- 
stream of the —10 element (5). This short sequence, the 



'extended —10' motif, is characterized by the 5'-TRTG-3' 
consensus sequence (6-9). Recognition of 'extended —10' 
promoters does not require specific contacts between 
region 4.2 and the —35 element (10,11). 

The galactose operon (ga/ETKM) of E. coli is trans- 
cribed from two overlapping promoters, PI and P2 
(Figure 1). PI is regulated by the galactose repressor 
(GaIR) by two different mechanisms, contact inhibition 
and DNA looping. GaIR dimers can bind to two operator 
elements, 0 E and 0\, separated by a 1 13 bp DNA segment 
comprising PI and P2 (Figure 1) (12-14). Contact inhib- 
ition occurs when 0 E -bound GaIR represses PI by con- 
tacting the C-terminal domain of the a subunit of RNAP 
(a-CTD), and inhibiting open complex formation (15-17). 
DNA looping repression of both PI and P2 occurs sim- 
ultaneously when the <9 E -bound GaIR dimer interacts 
with the Oj-bound GaIR dimer to form a tetramer and 
the intervening DNA loops out, forming a repression 
complex (known as repressosome) (18-20). Assembly of 
the Gal repressosome, a higher order nucleoprotein struc- 
ture containing an antiparallel DNA loop (21,22), requires 
(i) binding of two dimeric GaIR to 0 E and 0\, (ii) nega- 
tively supercoiled DNA, (iii) optimal angular orientation 
of 0 E and 0\ and (iv) specific binding of the architectural 
histone-like protein (HU) to a HU binding site (hhs) in the 
interoperator region (23,24). DNA supercoiling and the 
binding of HU at a critical position facilitate DNA loop 
formation over a short distance by decreasing the persist- 
ence length of DNA (25). 

Previous studies suggested that simultaneous binding of 
RNAP to the —10 and —35 elements introduces a specific 
bend in the DNA, as indicated by DNasel hypersensitivity 
around position —25 (4,26). In galPi, this hypersensitivity 
was less pronounced in the absence of a specific —35 rec- 
ognition site, suggesting that the DNA trajectory in 
RNAP-'extended —10' promoter complexes is different 
from that found in complexes where RNAP specifically 
recognizes both the —10 and —35 elements (11,26). 
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Figure 1. (Top) Schematic map of the gal regulatory region containing 
the gal operator sites (O e and Oi), the gal promoters, PI (+1) and P2 
(—5), and the HU binding site (hbs, +6.5). The transcription start point 
of PI is used as a reference in the numbering system (+1). (Bottom) 
V-shaped, stacked interaction of operator-bound GalR dimers can lead 
to two different antiparallel DNA loop trajectories, Al (left) and A2 
(right). Arrows indicate the direction of transcription initiated at the gal 
promoters. 



In this work, we asked whether the DNA bend 
introduced by simultaneous binding of RNAP to the 
— 10 and —35 elements affects the stability of the gal 
DNA loop. We introduced mutations at the —35 
hexamer of the PI promoter and monitored transcription 
regulation of the wild-type (WT) and mutant promoters in 
the absence and presence of GalR and HU by using 
in vitro transcription assays. Our results are consistent 
with a model where RNAP binding to the —10 and —35 
elements have an inhibitory effect on looping mediated 
repression of PI. 



MATERIALS AND METHODS 

Strain and plasmid constructions 

Plasmid manipulations followed protocols described in 
Sambrook and Russel (27). Transformations were per- 
formed with XL-1 Blue competent cells (Stratagene). 
Sequencing (ABI Prism) kits were purchased from 
Applied Biosystems, restriction endonucleases from Fer- 
mentas, DNA oligonucleotide primers from Invitrogen 
and DNA purification kits from Qiagen. DNA sequencing 
reactions were analysed in a Perkin-Elmer/Applied 
Biosystems (model 373A) automated sequencer. 

Mutations in the gal promoter regions were created 
by polymerase chain reaction (PCR) using the Platinum 
High Fidelity PCR SuperMix (Invitrogen) and inserted 
between the EcoRI and PstI sites in plasmid pSA850 
(28). The mutated sequences of the gal regulatory region 
in the resulting plasmids were verified. 

Protein purification 

Expression and purification of the hexahistidine-tagged 
GalR followed the protocol described by Semsey et al. 
(20). HU protein was purified according to the method 
described by Aki et al. (18). RNAP was purchased from 
USB. RNAP concentration was specified by the manufac- 
turer. GalR and HU concentrations were measured using 
the Micro BCA Protein Assay Kit (Pierce). The quality of 
protein preparations was tested in in vitro transcription 



reactions using the reference plasmid pSA850. Similar 
results were obtained to previously published results 
(22,23). 

In vitro transcription assays 

In vitro transcription reactions were performed as 
described in (20). The reaction mixture (50 ul) contained 
20 mM Tris acetate, pH 7.8, 10 mM magnesium acetate, 
200 mM potassium glutamate and 2 nM supercoiled DNA 
template. GalR concentrations vary from 5 to 40 nM as 
indicated, and HU was used at 80 nM. RNAP (20 nM) 
was added before incubating the reactions at 37°C for 
5min. Transcription was initiated by the addition of 
1.0 mM ATP, 0.1 mM GTP, 0.1 mM CTP, 0.01 mM 
UTP and 5uCi of [a- 32 P]UTP (3000Ci/mmol). 
Reactions were terminated after lOmin at 37°C by the 
addition of an equal volume of transcription loading 
buffer (0.025% bromophenol blue, 0.025% xylene 
cyanol, 0.01 M ethylenediaminetetraacetic acid and 90% 
deionized formamide). After heating at 90°C for 3 min, the 
samples were loaded onto an 8% polyacrylamide-urea 
DNA sequencing gel. RNA bands were quantified using 
the ImageQuant™ Phosphorlmager (Molecular 
Dynamics, CA). We followed the standard procedure 
that uses the RNA1 transcript as an internal control 
between lanes, to decrease the level of potential experi- 
mental error (19). The RNA1 transcript is not affected 
by GalR binding. Band intensities were background cor- 
rected as described previously (19). This procedure has 
<10% error (29). As levels of the studied transcripts 
relative to the level of the RNA1 transcript may slightly 
vary depending on the quality of the plasmid DNA prep- 
aration, promoter activities in the presence of GalR were 
expressed relative to the promoter activity in the absence 
of GalR. 

Construction of a mathematical model of the 

in vitro system 

In vitro transcription reactions contained a fix amount of 
GalR (0^10 nM), RNAP (20 nM) and DNA (2 nM). GalR 
has two specific binding sites (0 E and Oi) on the DNA 
molecules used, therefore, the amount of GalR added to 
the reaction (i? to tai) is the sum of the operator bound mol- 
ecules (R B ) and free GalR molecules (R F , not bound to 
DNA or bound non-specifically). As the DNA template is 
present at 2nM concentration, and there are two oper- 
ators per DNA molecule, i? B <4nM. Based on the 
possible patterns of operator occupancy, there are five 
different binding states for GalR: (i) not bound to any 
operators, (ii) bound to Oe but not to Oj, (hi) bound to 
Oi but not to 0 E , (iv) bound to both operators without 
DNA loop formation and (v) GalR dimers bound to 0 E 
and 0\ form a tetramer, resulting in a DNA loop. At any 
given GalR concentration, the DNA molecules present in 
the reaction are distributed between these five states. The 
ratios of these states can be computed based on the 
operator binding affinities and on the probability of 
DNA looping (Table 1). 

In each state, RNAP can specifically bind to only the 
— 10 region of the PI promoter, to both the —10 and —35 
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Table 1. The 15 binding states of the PI promoter and regulatory 
region considered in the model 



Stale 



0 E —35 —10 0i Loop Activity Relative weight 



1 

2 
3 
4 
5 
6 
7 
g 
9 
10 
11 
12 
13 
14 
15 



0 
I 
1 

0 

0.5 

0.67 

0 

I 

I 

0 

0.5 

0.67 

0 

0.05 
0.05 



Xf ;< Kf 
Kile Kri 



x >■ x L. 

Are Arj 

e- F x >, F x L x 

Are Ari *10 

Are * Ari * A'ioxF], * 



Parameters used are as follows: R F , concentration of free GalR; K RE , 
binding affinity of GalR to O e ; K m , binding affinity of GalR to Oi; L, 
looping factor; L , looping factor when RNAP is bound to both the 
— 10 and —35 elements; P, RNAP concentration; K Kh binding affinity 
of RNAP to the —10 region; Fj 5 , factor for increase in RNAP affinity 
to the promoter resulted from specific contacts at the —35 region. The 
five GalR-binding categories (i-v) are indicated on the left. 



regions, or be not bound to the PI promoter, resulting in 
15 binding states (Table 1). The relative statistical weights 
for the different binding states depend on the probabilities 
of GalR and RNAP binding to the operators and 
promoter elements, respectively. The weight of the state 
where DNA is not bound by any of the proteins is used as 
a reference (=1). The weights also depend on the prob- 
ability of DNA loop formation mediated by tetramer for- 
mation of GalR dimers bound to 0 E and 0\. Therefore, 
the relative statistical weights are functions of free protein 
concentrations and protein binding affinities for the dif- 
ferent sites on the DNA (Table 1). Under the quasi- 
steady-state assumption that GalR and RNAP binding 
to DNA are sufficiently rapid relative to the transcription 
rate, the probability of every 1 of the 1 5 binding states can 
be calculated (30). The probability of a state i (Pj) is equal 
to the weight of state i divided by the sum of the weights of 
all possible states. Promoter activity levels can be given by 
the following equation: 

Promoter activity = P, x activity,-. 

Relative promoter activities in each binding state (Table 1 ) 
are derived from experimental results. In the absence of 
RNAP binding, the promoter activity is 0, in the presence 
of DNA looping a weak promoter activity is observed 



(0.05), and contact inhibition by 0 E -bound GalR results 
in 50% repression of WT PI where RNAP does not make 
specific contacts with the —35 region (22) (relative activ- 
ity = 0.5), and 33% repression in the presence of the con- 
sensus —35 box (relative activity = 0.67). The latter value 
was determined by in vitro transcription assays performed 
on a pSA850 derivative plasmid containing the consensus 
PI -35 element (P1C35), in the presence of 80 nM GalR 
T322R, a GalR mutant that binds the gal operator sites 
similar to WT GalR but unable to form a DNA loop (19). 
Promoter activity of 1 was used in the RNAP-bound 
states where the 0 E operator is not bound by GalR. 
There are several other promoters on the plasmid 
template (28); hence, it is difficult to estimate the concen- 
tration of free RNAP. However, this problem can be 
overcome by using the term PjK w , where RNAP concen- 
tration is expressed relative to its binding affinity to WT 
PI. The upper limit for the value of P/K w is 0.02 because 
PI activity is increased 50-fold with the consensus —35 
box mutation (P1C35). We used this maximal value in 
the model, assuming that the P1C35 promoter is 
transcribed in all the template molecules in the in vitro 
transcription reactions. Concentration of R F was com- 
puted based on R to tai and R B , which were derived from 
the statistical weights listed in Table 1 . The mathematical 
model was written in FORTRAN. 



RESULTS AND DISCUSSION 

Effect of mutations in the —35 region on PI activity 

To study the role of the —35 element on PI regulation, a 
set of PI derivatives was created by PCR mutagenesis 
where the —35 PI element was randomized (YWNNNH, 
TNGNCA and TTGNNN). The mutant promoters were 
inserted into plasmid pSA850 (28), a reference plasmid 
used for in vitro studies on PI and P2 regulation. From 
the pool of plasmids containing mutant promoters, we 
selected 18 plasmids randomly and determined the 
sequence of the mutant promoter regions. The promoter 
set did not contain a promoter with a consensus —35 box; 
therefore, we created this PI derivative (TTGACA) sep- 
arately. Activity and regulation of the WT and 13 mutant 
promoters were studied using in vitro transcription assays 
(Figure 2). The same batches of RNAP, GalR and HU 
were used for all the assays to make the results compar- 
able. Mutations in the —35 box substantially increased the 
intrinsic activity of PI and decreased that of P2. The 
promoter containing the consensus —35 sequence 
showed the highest activity; it was ~50-fold stronger 
than WT PI (Figure 2, lane 22 versus lane 1; also refer 
Table 2). 

Typically the —35 elements of strong promoters are 
more similar to the consensus sequence than the —35 
elements of weaker promoters. However, less deviation 
from the consensus does not necessarily mean higher 
activity. For instance, a promoter with only two 
mismatches compared with the consensus sequence 
(TcGgCA) does not show substantially stronger activity 
than the WT promoter with six mismatches (cacttt) 
(Figure 2, lane 28 versus lane 1). 
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GaIR 
HU 

P2 
P1 



+ + 
- + 



+ + 
- + 



+ + 
- + 



+ + 
- + 



+ + 
- + 



+ + 
- + 



+ + 
- + 



RNA1 



1 2 3 
CACTTT 



4 5 6 
CTGGGA 



7 8 9 

TACGGA 



10 11 12 
TAGTGT 



13 14 15 

TTGTTA 



16 17 18 

TTGATG 



19 20 21 
TGGGCA 



GaIR 
HU 

P2 
P1 



RNA1 



+ + 
- + 



22 23 24 

TTGACA 



+ + 
- + 



25 26 27 
TTGTGT 



+ + 
- + 



28 29 30 

TCGGCA 



+ + 



+ + 



+ + 



+ + 




40 41 42 
TGGACA 



Figure 2. Results of in vitro transcription assays performed on WT (lanes 1-3) and mutant —35 derivatives (lanes 4-42) of PL GaIR and HU 
concentrations were 40 and 80 nM, respectively. The RNA1 transcript, which is not affected by GaIR and HU binding, serves as an internal control 
between lanes. 



Table 2. Relative activities of WT and mutant promoters in the 
presence of GaIR and 80 nM HU 



Promoter Strength Relative activity 

(P/RNA1) GaIR 







0 nM 5 nM 


10 nM 


20 nM 


40 nM 


PI 


0.2 


1 0.27 


0.09 


0.08 


0.06 


P7C35 


10.0 


1 0.71 


0.36 


0.22 


0.15 


WS35G 


4.1 


1 0.72 


0.38 


0.16 


0.06 


WS35C 


2.6 


1 0.74 


0.42 


0.20 


0.07 


WOO 


2.5 


1 0.52 


0.18 


0.09 


0.06 


P2 


0.3 


1 0.54 


0.17 


0.09 


0.07 


P2C35 


1.9 


1 0.89 


0.50 


0.18 


0.05 



Promoter strength is expressed as the ratio of intensities of bands cor- 
responding to the unregulated promoter and the RNA1 transcript. 
Promoter activities were determined based on the results of in vitro 
transcription assays presented in Figure 3, as described in the 
Materials and Methods section. Concentrations are calculated for 
GaIR dimers. 



Effect of mutations in the —35 region on the 
regulation of PI 

Mutations in the —35 box of PI did not result in a quali- 
tative change in the regulation of promoters. Similar to the 
WT promoters, mutant PI promoters were repressed in 
the presence of GaIR, while P2 promoters were activated. 
In the presence of both GaIR and HU, both the PI and 
the P2 promoters were repressed as a result of DNA loop 
formation (Figure 2). The results suggested that the effi- 
ciency of PI repression depends on the —35 sequence. 
HU-independent GalR-mediated transcription inhibition, 
which is a combined result of contact inhibition and 
non-assisted DNA looping (22,31), was least efficient in 
the case of the promoter containing the consensus (TTGA 
CA) -35 sequence. This promoter retained 50% activity 



in the presence of GaIR, compared with 20% activity 
observed in the case of the WT promoter (Figure 2, lane 
23 versus lane 22 and lane 2 versus lane 1). Promoters with 
TgGgCA and TTGttg —35 hexamers showed similar 
level of GalR-mediated repression as the WT promoter 
(Figure 2, lanes 19-20, 34^35). Results obtained in the 
concomitant presence of GaIR and HU suggested that 
Gal repressosome-mediated repression can also be 
affected by certain mutations in the —35 boxes. Repres- 
sion of most of the mutant promoters by repressosome- 
mediated DNA looping was comparable with the WT 
promoter. However, in some cases we observed more in- 
complete repression, which was most prominent in the 
case of the promoter containing the consensus —35 
sequence (TTGACA) (Figure 2, lane 24 versus lane 22). 
Two other promoters, which had only one mismatch in the 
—35 hexamer compared with the consensus sequence 
(TgGACA and TcGACA), also showed incomplete 
repression (lane 42 versus lane 40, and lane 33 versus 
lane 31). 

Looping repression inhibits a step prior the formation 
of the first phosphodiester bond (32). Torsional and 
lateral inflexibility of the short looped DNA may impede 
RNAP binding or isomerization from close to open 
complex (33). The lifetime of the gal DNA loop is rela- 
tively short and decreases with increased tension in the 
DNA. The time spent in the unlooped state is also 
tension dependent and increases with tension (34). To 
obtain efficient repression, the time interval when the 
DNA is in the unlooped state must be shorter than the 
time RNAP needs to find the promoter (35) and initiate 
transcription. Therefore, the DNA loop should form 
faster if DNA looping inhibits only open complex forma- 
tion and not RNAP binding to the DNA. 

The length of the time interval spent in the unlooped 
state is not known in our in vitro transcription assays but 
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perhaps RNAP is allowed to bind to the gal promoters in 
the unlooped intervals. This assumption is based on the 
observation that the A2 loop trajectory (Figure 1 , bottom 
right), which has similar thermodynamic properties as the 
Al loop (formed on WT DNA), fails to repress transcrip- 
tion, likely because the RNAP binding site overlaps with 
the HU binding site in the A2 loop (22,34). In this case, 
RNAP binding inhibits reclosing of the loop with the A2 
trajectory. As the Al and A2 trajectories have similar life- 
times of looped and unlooped states (the A2 loop forms 
and breaks down slightly more frequently) (34), we can 
assume that RNAP can bind to the gal promoters in the 
unlooped intervals of the WT setup used in the in vitro 
transcription assays. 

There are at least three different possible explanations 
for the incomplete looping repression: (i) mutations alter 
the deformability of DNA; (ii) stronger promoters need 
shorter times for initiation of transcription and (iii) the 
bending introduced by concomitant binding of RNAP to 
the —35 and —10 elements interferes with DNA loop for- 
mation. In regard to (i), it is unlikely that the altered —35 
sequence would substantially change DNA deformability. 
The predicted twisting difference (36) caused by the TTG 
ACA substitution is very minor (~0.2°), and the substitu- 
tion is located in a region which is not deformed in the 
DNA loop (22). Therefore, in the next part of the manu- 
script, we focus on the other two scenarios. To study 
the second possibility (ii), we tested whether increasing 
the promoter strength by enhancing the —10 element has 
similar effect as observed in the case of —35 element mu- 
tations. The third possibility (iii) was tested by a compu- 
tational approach, asking whether the data obtained on 
the WT and on the -35 mutant (TTG AC A) PI promoters 
can be fitted to the model using the same GalR binding 
and looping parameters. 



The effect of promoter strength on looping repression 

As described earlier, looping repression becomes ineffi- 
cient if the unlooped time is shorter than the time 
required for transcription initiation at the promoter 
located inside the DNA loop. To test whether incomplete 
repression obtained for some of the mutant promoters 
(Figure 2) can be explained by the increased promoter 
strength (faster initiation), we created a mutant, which 
has a consensus —10 hexamer (Figure 3, P1C10). Repres- 
sion of the P1C10 mutant was compared with repression 
of the —35 hexamer mutants, which showed incomplete 
looping repression, TTGACA (P1C35), TgGACA 
(P1S35G) and TcGACA (P1S35C) at different GalR con- 
centrations (Figure 3 and Table 2). We did not find a 
strong correlation between promoter strength and effi- 
ciency of repression. The P1C10 promoter has similar in- 
trinsic activity to the P1S35C promoter; however, 
repression of P1S35C is less efficient at intermediate 
GalR concentrations. Also, while the intrinsic activity of 
P1S35G is somewhat stronger than that of P1S35C, re- 
pression of these promoters showed similar dependence 
on GalR concentration. A mutant P2 promoter having a 
consensus —35 box was also weaker then PI CIO but was 
less efficiently repressed by GalR than P1C10 or the WT 
P2 promoter (Table 2). Therefore, we concluded that 
faster transcription initiation by itself is not sufficient to 
explain the behavior of the —35 mutant promoters. 

An important difference between the P1C10 and 
the —35 mutant PI promoters is that RNAP binding to 
the —35 mutant promoters results in a specific bend in the 
DNA. Stability of short DNA loops is often increased by 
architectural proteins which bind to the looped region and 
introduce a bend with a characteristic stereo specificity 
(25). However, introduction of a bend at an improper 
position can inhibit DNA loop formation (37). The bend 



A WT (PI) tGTGTAAACGATTCCACtaatttattccatqtcacacttttcqcatctttcft TATG c TATqqT tatttcAta 

WT (P2) tGTGTAAACGATTCCACtaatttattccatqtcacacttttcqcatc TtTG t TATqcT atqqttAtttcata 

P1C35 tGTGTAAACGATTCCACtaatttattccatqtcaTTGACAtcqcatctttqt TATG c TATqqT tatttcAta 

P1S35G tGTGTAAACGATTCCACtaatttattccatgtcaTgGACAtcgcatctttgt TATG c TATqgT tatttcAta 

P1S35C tGTGTAAACGATTCCACtaatttattccatqtcaTcGACAtcqcatctttqt TATG c TATqqT tatttcAta 

P1C10 tGTGTAAACGATTCCACtaatttattccatgtcacacttttcgcatctttgtTATCcTATAATtatttcAta 

P2C35 tGTGTAAACGATTCCACtaatttattccaTTGACAacttttcgcatc TtTG t TATqcT atggttAtttcata 

0 E -35 -10 tsp 

B 

GalR (nM) 0 5 10 20 40 0 5 10 20 40 0 5 10 20 40 0 5 10 20 40 0 5 10 20 40 0 5 10 20 40 




WT P1C35 P1S35G P1S35C P1C10 P2C35 



Figure 3. (A) DNA sequences of PI, P2 and their mutant derivatives. The 0 E operator sequence (blue) was used to align the sequences. 
Transcription start points (tsp) are typed boldface. Extended —10 and —10 promoter elements are underlined. Red sequences indicate mutations. 
Positions matching the consensus —10 and —35 sequences are capitalized. (B) Results of in vitro transcription assays performed on plasmids carrying 
the indicated promoter sequences in the presence of O^tOnM GalR. HU (80 nM) was present in all of the reactions. The same batches of RNAP, 
GalR and HU were used for all the assays to make the results comparable. These batches were different from the ones used in the experiments shown 
in Figure 2. Band intensities were quantified and are shown in Table 2. 
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introduced by RNAP in the case of P1C35 (1 1) most likely 
interferes with DNA loop formation because it affects a 
region which is not deformed in the GalR-mediated DNA 
loop formed on WT DNA. Inhibition of DNA loop for- 
mation allows longer time for RNAP binding and tran- 
scription initiation, resulting in failure of DNA loop 
mediated repression as it was observed in the A2 loop 
geometry, where RNAP binding to the P2 promoter 
inhibits HU binding and stabilization of the DNA loop 
(22,34). 

Mathematical model of the in vitro system 

We have constructed a mathematical model to simulate 
the effects of specific RNAP contacts in the —35 region. 
The model is described in the 'Materials and Methods' 
section. By fitting this model to the experimental data 
on WT PI and P2 (Table 2), we obtained estimates on 
the parameters for the binding strengths of GalR to 0 E 
C^re) and 0\ (^ RI ), and also for the looping factor (L). 
The fitting criterion used was that the predicted values 
must not be >20% different from the corresponding 
measured values (Figure 4). This interval is at least twice 
as big as the expected experimental error (29). -F35 was set 
to 10000, assuming negligible contribution of specific 
RNAP binding to the —35 region. Parameter sets Crv RE , 
K K i and L) were chosen randomly and the sets fulfilling 
the fitting criterion were recorded. In the parameter sets 
obtained, ^T RE was between 1.1 and 7.6 nM. Previous es- 
timates of A^ RE by DNasel footprinting (1.3 nM) (38) and 
by fluorescence anisotropy (4.2 nM) (39) fall in this 
interval. The smallest values obtained for K m and L 
were 36 nM and 137, respectively. Based on these values, 
we can conclude that the states where GalR is bound to 0\ 
but not to 0 E or bound to both operators without DNA 
loop formation have little importance because 
Pf/K re » R F /K m , and L/K m » 1/K m . In the rest 
of the GalR binding states, the parameters L and K RE 

appear only in the looped state as k ^**k ■ Therefore, in 

the parameter sets that satisfy the fitting criteria the L/K m 
ratio is constrained, falling between 2.3 and 4.1 nM -1 . 

To fit the model to the data on the P1C35 promoter, we 
used F 35 = 0.0001 (assuming concomitant —10 and —35 
contacts), ^T RE values from the interval obtained on the 
WT promoters (1.1-7.6 nM), and _K RI >36nM. We again 
accepted parameter sets that give activities within 20% 
from the measured values (Figure 4). In the parameter 
sets obtained, AT RI values were between 36 and 500 nM, 
L values were between 1 and 136, whereas L/K Ki ratios 
were between 0.021 and 0.271 nM -1 . Data could be fitted 
using Fj 5 values that were smaller than 0.001. However, 
L values were always < 136, and L/K^i ratios were 
<0.52nM _1 . Therefore, it was not possible to fit the 
data obtained on the WT and on the —35 mutant (TTG 
ACA) PI promoters using the same GalR binding and 
looping parameters. Results of the simulations indicate 
that DNA loop formation is less probable in the case of 
the P1C35 promoter, i.e. DNA loop formation is in- 
hibited. Unlike the P1C35 data, results obtained on the 
P1C10 promoter could be explained by weaker repression 
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Figure 4. Fitting of the mathematical model to the experimental data 
obtained on WT PI (top), P2 (middle) and P1C35 (bottom). Crosses 
indicate experimental measurements (shown also in Table 2) and 
vertical bars show the ±20% acceptance intervals. Gray areas 
indicate results of simulations. i? to tai indicates the concentration of 
GalR dimers (sum of the free and operator bound forms). 



by the 0 E -bound GalR dimer. These data could be fitted 
using L/K K i >2 and 33% repression by contact inhibition. 
Therefore, it is likely that specific RNAP binding to the 
—35 sequence is required for inhibition of DNA loop 
formation. 

Using the parameters obtained, we computed the distri- 
bution of the different binding states as a function of total 
GalR concentration in the reactions. The left panel of 
Figure 5 shows how the DNA molecules carrying the 
WT gal regulatory region are distributed between the dif- 
ferent GalR-bound states. At low GalR concentrations 
the unbound state dominates, whereas at higher GalR 
concentrations, most of the DNA molecules are in the 
looped state. The concentration of single operator 
bound states shows non-monotonic behavior because of 
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Figure 5. Modeling the distribution of DNA molecules containing the WT PI promoter (left) and the P1C35 promoter (right) between the five 
different binding states of GalR as a function of GalR concentration (R Mll i). Hatched areas between the boundary lines indicate the concentration 
ranges of DNA molecules present in the different binding states. The five binding states of GalR are the following: (i) not bound to any operators 
(solid green boundary lines), (ii) bound to O e but not to O t (solid blue boundary lines), (iii) bound to 0\ but not to O e (dashed green boundary 
lines), (iv) bound to both operators without DNA loop formation (dashed blue boundary lines) and (v) bound to both operators and forming a DNA 
loop (solid red boundary lines). Same colors mean same regulatory outcome: no repression (green), contact inhibition by the 0 E -bound GalR dimer 
(blue), and repression by DNA looping (red). Vertical dotted lines at 1.5 and 15nM GalR concentrations are shown for guidance. In the schematic 
drawings of the binding states, (top) gray boxes indicate GalR, and the arrow indicates the direction of transcription initiated at the promoter 
studied. 



the cooperative interactions in DNA loop formation. The 
state where GalR binds to only 0\ is negligible, whereas 
~5-20% of the DNA molecules are in the only 0 E -bound 
state at 3nM GalR concentration (Figure 5). In the 
absence of D-galactose, the steady-state level of intracellu- 
lar free GalR (R F ) is ~15nM (40). At the corresponding 
GalR concentration (i? t otai — 19nM) >95% of the mol- 
ecules are in the looped state. 

We performed a similar analysis of states in the case of 
the P1C35 promoter. The curves had similar trajectories to 
the corresponding ones obtained with the WT promoter, 
however, they were quantitatively different (Figure 5, right 
panel). At R F = 15 nM (i? total -18.5 nM) only -75% of 
the molecules were in the looped state, and —10-20% of 
molecules were in the state where only 0 E is bound. The 
highest abundance of the latter state (15-50%) was 
observed at i? to tai of 4-5 nM. 

The weaker loop formation in the presence of specific 
RNAP contacts with the —35 region may have important 
consequences in vivo. The gal DNA loop is not able to 
repress transcription after elongation is initiated. 
Therefore, a single elongation initiation results in dissoci- 
ation of the Oi-bound GalR dimer and therefore increases 
the time spent in the unlooped state. This event can facili- 
tate multiple initiations because the on-rate of intracellu- 
lar free RNAP (present at few hundred nanomolarity 
concentration) (41) is much higher than the on-rate for 
GalR binding to 0\ and then forming a loop. This effect 
can become more pronounced when D-galactose is present 
at low levels, decreasing the active fraction of free GalR. 
Besides inhibiting DNA binding, the presence of D-galact- 
ose also interferes with GalR tetramerization resulting in 
an increase of the fraction of time spent in the unlooped 
state (42). Interaction of GalR dimers is more sensitive to 



D-galactose concentration than operator binding by GalR. 
For instance, at 0.1 mM D-galactose GalR tetramerization 
is completely inhibited, whereas —85% of the 0 B and 
>50% of the Oi operator DNA can be bound by GalR 
(42). Therefore, burstiness of transcription initiation 
resulted from the slow on-rate of a transcription regula- 
tory protein (35) would likely modulate the behavior of 
the feedbacks circuits in the galactose utilization system 
(43,44) in a D-galactose-dependent manner. There are two 
mechanisms that act against such burstiness of transcrip- 
tion initiation at the galPl promoter: lack of specific 
RNAP contacts at the —35 region, and inhibition of 
open complex formation by the 0 E -bound GalR dimer 
(17). The first mechanism reduces the interference of 
RNAP binding and DNA loop formation, whereas the 
second mechanism decreases the probability of transcrip- 
tion initiation, and thus the chance of removal of GalR 
from 0\ by the elongating RNAP. 



CONCLUSIONS 

The extended —10 motif is present in —20% of all E. coli 
promoters (6), and 26% of these promoters have weak 
—35 boxes, with not more than two matches to the con- 
sensus sequence (8). Promoters belonging to this minor 
class of promoters lacking the specific —35 element recog- 
nition are utilized to transcribe the gal operon. RNAP 
binding to these promoters does not interfere with forma- 
tion of the short gal DNA loop, suggesting that it has little 
if any impact on DNA flexibility. However, further studies 
are needed to investigate whether extended —10 promoters 
were generally evolved to decouple local chromosome 
structure from transcription initiation. 
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